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ABSTRACT 

We study how X-rays from stellar binary systems and the hot intracluster medium 
(ICM) affect the radiative cooling rates of gas in galaxies. Our study uses a novel 
implementation of gas cooling in the moving-mesh hydrodynamics code AREPO. X- 
rays from stellar binaries do not affect cooling at all as their emission spectrum is too 
hard to effectively couple with galactic gas. In contrast, X-rays from the ICM couple 
well with gas in the temperature range 10'^ —10® K. Idealised simulations show that the 
hot halo radiation field has minimal impact on the dynamics of cooling flows in clusters 
because of the high virial temperature (> lO^K), making the interaction between the 
gas and incident photons very ineffective. Satellite galaxies in cluster environments, 
on the other hand, experience a high radiation flux due to the emission from the 
host halo. Low mass satellites (< IO^^Mq) in particular have virial temperatures 
that are exactly in the regime where the effect of the radiation field is maximal. 
Idealised simulations of satellite galaxies including only the effect of host halo radiation 
(no ram pressure stripping or tidal effects) fields show a drastic reduction in the 
amount of cool gas formed (~ 40%) on a short timescale of about 0.5 Gyrs. A galaxy 
merger simulation including all the other environmental quenching mechanisms, shows 
about 20% reduction in the stellar mass of the satellite and about ~ 30% reduction 
in star formation rate after 1 Gyr due to the host hot halo radiation field. These 
results indicate that the hot halo radiation fields potentially play an important role 
in quenching galaxies in cluster environments. 

Key words: atomic processes - hydrodynamics - plasmas - radiative transfer - 
galaxies: formation - methods: numerical 


1 INTRODUCTION 


The basic model of galaxy formation assumes that as dark 
matter hierarchically collapses, gas flows into gravitational 
potential wells and forms stars at their centers ( [White ^ 
Rees|1978 1. The gas converts gravitational potential energy 


into thermal energy and consequently, the gas in the halo 
heats up to the virial temperature, r„ir. In the most massive 
galaxies, the temperature becomes so large that the cooling 


time for the gas is longer than the age of the Universe (Rees 
jfc Ostriker||1977| ). In this idealised picture, gas in massive 
galaxies cannot cool to form stars. However, galaxy haloes 
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have high density material in their centers, and since gas 
cooling is strongly dependent on density, A oc n^, this dense 
central gas unstably overcools and collapses to the halo cen¬ 


ter where it forms many more stars than are observed (Cole 
1991 Navarro et al.||1995 Benson et al.|2003 l. 


To counteract this overcooling, galaxy formation mod- 


els feed energy back from stars (e.g.. 

Navarro et al. 19961 

Springel & Hernquist|2003||Stinson et a 

.[2006 Dalla Vecchia 

& Schaye 2008 Vogelsberger et al. 2013 

Hopkins et al. 

20141 

and active galactic nuclei (AGN) (e.g., [Springel et al. 

2005 


Di Matteo et alT]|2005| jSijacki et al.||^07[ |2009| [Booth fcj 


Schaye 20091. Stars deposit large amounts of energy through 


stellar winds, supernovae explosions, and radiation (Lei- 
therer et al.[[1999l. Although feedback from stars is able to 
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2 Kannan et al. 


regulate the star formation rate in low mass galaxies, they 
are inefficient at high masses ( Kannan et al.|[2014a l. AGN 
feedback has been invoked to solve the problem in high mass 
haloes (Mhaio ^ 10^^M*).The details of how AGN couple to 
the gas in and around galaxies is still uncertain, so modeling 
efforts have so far been necessarily crude. 

Springel et al. ( 2005|), for example, used the energy from 


need to be independent of the star formation in galaxies. 
Either the stellar emission sources must have a time delay 
between the star formation event and emission or they must 
be of a non-stellar origin. 


Gnedin & Hollon (20121 created a general model for 


the high accretion periods of AGN when they shine as quasi- 
stellar objects (QSOs, quasars) to heat gas immediately sur¬ 
rounding the black hole and drive large scale winds. How¬ 
ever, this level of accretion requires that a lot of cold gas 
is present at the center of the galaxy, which in turn leads 
to star formation. Even when AGNs are accreting a small 
amount of gas, they drive spectacular radio jets that are 
well collimated by magnetic fields and produce radio lobes 
tens of kpc away from the black holes. Gosmological numer¬ 
ical simulations have cited these jets as a form of feedback 
that is able to heat the entire gas halo to high enough tem¬ 
peratures that prevents gas from cooling into the disk and 


cooling in the presence of a radiation field near a galaxy 
(including both stars and AGNs). They showed that these 
radiation fields will significantly alter the cooling rates of gas 
in the vicinity of a massive ’O’ star and in galaxies contain¬ 
ing highly luminous quasars. Along these lines |Vogelsberger| 
et al. (20131 introduced a AGN feedback mechanism, which 


forming stars (Sijacki et al. 20071. However, these models 


do not address one of the primary features of radio mode 
feedback, the strong directionality of the jet. So neither of 
these models sufficiently explain how galaxies are quenched 
by AGN feedback. 

Thermal conduction is another mechanism which has 


modified the cooling of gas according to the incident ra¬ 
diation field of the central AGN. In their simulations, the 
radiative feedback is only effective during times when the 
AGN is strongly accreting and shining brightly as a quasar, 
which again corresponds to times when stars are forming. In 
the end, the only feedback that they found effective at main¬ 
taining quiescent galaxies was the AGN radio-mode feedback 
that operates by heating the halo gas to temperatures where 
its cooling time is too long to cool into the disk to form stars. 

Another avenue through which galaxies can quench is 
through environmental effects such as ram pressure strip¬ 
ping ( Gunn fc Gott||1972 McGarthy et al.|2008 Font et al. 


2008 Kang fc van den Bosch||2008 Weinmann et al.||2010 1 
and tidal stripping (Taylor fc Babul||2001 D’Onghia et al. 


N arayani 

2003 

IVoit et al. 1120081 [20141 |Voit||2011| |Voit & 

invoked to explain the different statistical properties of field 

Donahue 

2015 

1. The large Ghandra study of ICM entropy 

and cluster galaxies, such as the prevalence of ellipticals and 


profiles (ACGEPT), showed that strong AGN feedback is 


present only in clusters with low-entropy cores (Cavagnolo 
et al.|200^ 20091. Voit et al. (2014 \ have postulated that the 


low entropy clusters have cooling times about the order of 
precipitation times of the clusters, which is empirically seen 
to occur when tcooi ~ lOtg (Gaspari et al.|2011 20121. The 


lower star formation rates and hence redder colors. In many 
prescriptions for ram pressure stripping the hot gas in the 
halo is assumed to be stripped immediately after the satel¬ 
lite falls into the host, which leads to a rapid decline in star 
formation rates (Baldry et al.|[2006 Weinmann et al!||2006 


clusters with large entropy at their centers were shown to 


Wang et al. 2007). More improved methods of stripping have 
been looked at, such as the a more gradual stripping of the 


lie above the thermal conduction balance solution, meaning 

hot gas reservoir (Font et al. 

2008 Kang & van den Bosch 

that the cooling in the centers of these clusters is balanced 

2008 Weinmann et al.||2010 

McCarthy et al. 20081. How- 


by thermal conduction heating by hot gas at larger radii. 
One of the most important assumptions in this study is the 
value of the conduction coefficient which is about 30% the 
canonical Spitzer conductivity ( Spitzer|1962|). Thermal con- 


ever, even this formalism does not agree well with the ob¬ 
servational measurements and they generally overproduce 
the quenched fraction of satellites as a function of distance 
from the cluster and underproduce the fraction of elliptical 


ductivity can reach this high value only if the magnetic field 
is chaotic over a wide range of length scales (factor of 100 or 


more), as might happen with MHD turbulence (Narayan & 


galaxies in cluster environments (Guo et al. 2011 Weinmann 
et al.|2011|). 


Medvedev 20011. However, recent studies have shown that 


the reproduction of observed sharp features in the ICM (cold 
filaments, bubbles) require that the conduction coefficie nt be 
as low as 10“^ — 10“^ times the Spitzer conductivity (Gas- 
W3l|Arth et al.|2014|). 


pari & Churazov 


In this paper we explore the possibility that local X-ray 
radiation from X-ray binaries (XRBs) and bremsstrahlung 
cooling radiation from hot gas present in the haloes of galax¬ 
ies can help to regulate star formation in galaxies. XRBs are 
known to be the most prominent X-ray sources in a galaxy 
which is devoid of an accreting black hole or large amounts 


Cantalupo (20101 explored analytically the effect of of hot optically thin gas (Fabbiano fc White 20031. Glus- 


X-ray radiation from hot shock heated gas due to stellar 
winds and supernovae (SNe) explosions on the cooling of 
halo gas. The high-energy radiation from these sources was 
shown to ionize metals and decrease the cooling rate of high- 


ters of galaxies are the brightest extended X-ray sources in 

the Uni verse, with their luminosities ranging from 

erg s“^ (Vikhlinin et al.||2009 Pratt et al. 2009 Anderson 


metallicity gas considerably. Kannan et al. (2014b hereafter 


RK14b) implemented this effect using GASOLINE (Wadsley 
et al.|2004 1 and showed that in a cosmological simulation of 


a Milky-Way galaxy the local radiation fields indeed reduce 
gas cooling rate onto the disc of the galaxy thereby reducing 
the star formation. This formalism will not quench galaxies 
because of the need to form stars in order to reduce star 
formation. The radiation sources which can quench galaxies 


et al. [2014 1. We build on the RK14b model and account for 

these additional X-ray radiation fields while calculating gas 
cooling rates in galaxy formation simulations. 

The paper is structured as follows. Section describes 
the photoionization sources that we consider in this work. 
Section outlines the approximations used in our radiative 
transfer approach, while Section|^describes the construction 
of the cooling table. Finally, in Section and we present 
the results of our implementation of the local photoioniza- 
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Local Photoionization Feedback -II 3 



Figure 1. The incident photon flux as a function of energy, from 
high (green curve) and low (cyan curve) mass X-ray binaries com¬ 
pared to the radiation field from star formation events or young 
stars (blue curve) and from old or post-AGB stars (red curve). 
These fluxes are reported at 10 kpc from the source. The radia¬ 
tion flux from both the HMXBs and the young stars are quoted 
for a star formation rate of 10 Mq yr“^. The radiation flux 
from both old stars and LMXBs are quoted for a stellar mass 
of M* = 1O^^M0. For comparison, the estimate for the flux of 
the extragalactic UV background at 2 = 0 is also plotted (purple 
curve). 


tion feedback (LPF) model on a test gas particle and on 
isolated galaxy simulations. Our conclusions are presented 
in Section 0 


2 IONIZING SOURCES 


There are several radiation sources that produce high energy 
ionizing photons. RK14b modelled the effect of radiation 
from star formation events (or young stars) and post-AGB 
(old stars) stars. In this paper we explore the emission from 
X-ray binaries (XRBs) and the cooling radiation from the 
hot tenuous gas around haloes. Throughout this work, the 
estimate for the extragalactic UV background by [Faucher^ 


Giguere et al. (20091 is used. 


2.1 X-Ray Binaries 


About a third of all stars form as part of a binary system 
I Lada|[^006 1 . The evolutionary timescales of stars compris¬ 
ing the binary system may be different, leading to the col¬ 
lapse of one star earlier than the other. The collapsed star 
then acts as an accretor of material from the donor main 
sequence star, once the surface of the donor star crosses 


the inner Lagrange point of the system (Paczyhski 19711. 


This infalling material dissipates gravitational potential en¬ 
ergy in the form of radiation. They are the most prominent 
X-ray sources in galaxies which are devoid of an accreting 


black hole or large amounts of hot optically thin gas ( Fab- 
biano &: White|[2003 1. The XRB population is divided into 
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Figure 2. The volumetric cooling (solid curves) and heating (dot 
dashed curves) rates of gas under the influence of radiation from 
young stars and HMXBs (blue curves), only young stars (green 
curves), old stars and LMXBs (cyan curves), only old stars (red 
curves) and the extragalactic UV background (purple curves). 
The gas is at a density of nu = 0.001cm“®, and a metallicity of 
Zq. The radiation flux for each component is equal to the value 
quoted in Fig. 


three categories depending on the mass of the donor star: 
high mass XRBs (M* is 8 M©), intermediate mass XRBs 
(1 Mq < M* < 8 Mq) and low mass XRBs (M* ^ 1 Mq). 
The X-ray emission is dominated by high mass XRBs im¬ 
mediately after a star formation event, followed by a possi¬ 
ble contribution from the intermediate mass XRB popula¬ 
tion, and the late time emission is dominated by low mass 
XRBs. The emission from and properties of, intermediate 
mass XRBs is not well studied because there are very few 


XRBs in this mass range in the Milky-Way galaxy (Podsi- 


adlowski et al. 2002 Pfahl et al. 20031. Therefore, in this 


work we only consider the emission from the high and low 
mass XRB populations only. 


2.1.1 High Mass X-ray Binaries 


High mass XRBs (HMXBs) are binaries where the mass of 
the donor star is > 8 M 0 . Their lifetimes are determined by 
the nuclear timescale of the high mass donor star, which is 
1 — 10 Myr. Thus, HMXB emission occurs a few million years 
after star formation events and the total emission luminosity 
of HMXBs follows the star formation rate of the galaxy. 
Tight empirical relations are observed between the X-ray 
luminosity of the HMXBs and the star formation rate of 
the galaxy ( Grimm et al.|20'03 Persic et al. [2004 1. Recently 
Mineo et al. (20121 presented the following relation for the 


luminosity of the HMXBs as a function of the star formation 
rate 


= 2.61 X 10^® 


SFR 


(iM0yr-i 


erg s 


( 1 ) 


Energies lower than 0.5 keV are difhcult to observe, because 
they are absorbed in the interstellar medium of the host 






































4 Kannan et al. 


as well as our own Milky-Way galaxy. We therefore model 
XRB emission in the hard X-ray regime only. In this regime 
(above 0.5 keV), the HMXB SED can be approximated by 
a power law: 


dl 

dE 


oc 


( 2 ) 


where I is the intensity and E is the energy of the emit¬ 
ted photon. Using these relations we can plot the SED of 
HMXBs as a function of the star formation rate of the 
galaxy. Fig. a shows the HMXB emission spectra corre¬ 
sponding to a SFR of 10 Mq yr“^ (green curve). For com¬ 
parison, the SED from young stars (see RK14b for details) 
is shown for the same SFR (blue curve). The fluxes quoted 
are at a distance of 10 kpc from the source. 

Figj^ shows the corresponding cooling (solid curves) 
and heating (dashed curves) rates in presence of various 
radiation sources. The gas parcel is at a density of tih = 
0.001cm“^ and Z = Z© and at a distance 10 kpc from the 
source. The cooling curve is calculated under the assump¬ 
tion that the gas in between the source and the gas par¬ 
cel is optically thin. Thus, the flux is inversely proportional 
to the distance squared. The heating (H) and cooling (A) 
curves of the gas in the presence of the young stars (blue 
curves) radiation field only lies exactly on top of the heating 
and cooling curves calculated in the presence of young stars 
-I- HMXB radiation field (green curves). In a star forming 
galaxy the emission from young stars (i.e., blackbody emis¬ 
sion from stars and X-ray emission from shock heated gas) 
dominates the emission from HMXBs up to ~ 180 Rydbergs 
(Fig.g. Above this energy, the hard X-ray photons do not 
couple effectively with the gas because of low interaction 
cross section. We conclude that HMXBs, although an im¬ 
portant source of X-rays have a sub-dominant effect on the 
gas cooling rates in galaxies. 


2.1.2 Low Mass X-Ray Binaries 


Low Mass XRBs (LMXBs) are binaries where the mass of 
the donor star is < IM©. The X-ray active phase of these 
systems are delayed by the nuclear timescale of the low 
mass donor star or by the orbital decay timescale. These 
timescales are ~ 1 — 10 Gyrs (Verbunt & van den Heuvel 
|1995[ ). The duration of the X-ray active phase also has simi¬ 
lar timescales | Podsiadlowski et al.|200'2 Pfahl et al.|2003 |. 
This delay in X-ray active phase decouples the X-ray emis¬ 
sion from the star formation event. Instead, the luminosity 
of LMXB correlates well with the total stellar mass of the 
galaxy ( Gilfano^|2004 1, and the relation is given by: 


= 1.00 X 10"° (iS^) 

The spectral shape of LMXB emission is assumed to 
follow the same functional form as Eq. i Fig. 0 shows the 
LMXB emission flux (cyan curve) for a galaxy with a total 
stellar mass of IO^^Mq at a distance of 10 kpc from its 
center. For comparison the flux from post-AGB stars (red 
curve), under similar conditions is also shown. The LMXBs 
contribution is primarily in the hard X-ray regime, while 
post-AGB stars shine only in the UV. 

The corresponding cooling curves are plotted in Figj^ 


The cooling curves for a parcel of gas under the influence of 
radiation from just post-AGB stars (red curve) and one with 
radiation from post-AGB and LMXBs (cyan curve) are in¬ 
distinguishable, except for a slight suppression of cooling at 
log^Q(T[K]) ~ 5.3 in the presence of LMXB radiation field. 
This implies that XRBs in general do not have a significant 
impact on the cooling rates. This is due to a couple of fac¬ 
tors, one, their emission fluxes are sub dominant compared 
to the radiation sources already considered in RK14b, i.e., 
young stars and post-AGB stars. And two, the hard X-ray 
photons have very low interaction cross sections. If we could 
in fact include the XRB emission at lower energies, then it 
could possibly have an impact on gas cooling rates. How¬ 
ever, modelling this part of the emission spectra has so far 
proven to be quite challenging (Gilfanov fc Bogdan||2010 1. 


2.2 Hot Gas Emission 


Early observations indicated that clusters are ubiqui¬ 
tous emitters i n the X-ray regime , with luminosities ~ 
104°“45g]-g g-i ( Giacconi et al.|l972 1. The X-ray emission is 
observed to be non time-varying and extended, with the spa¬ 


tial extent matching the galaxy distribution (Kellogg et al. 


1972 Elvis 19761. These properties point to the fact that 


the X-ray emission arises from bremsstrahlung or free-free 
cooling of the hot (T~ 10^“®K) tenuous (uh ~ 10“° cm“°) 
intracluster medium (ICM). The thermal free-free volumet¬ 
ric emissivity {jl^) of a parcel of gas at temperature T is 
given by 

— A^{T)neni, 


J 

where 


A„(r) = 6.8 X 


(4) 


(5) 


where p// is the free-free gaunt factor, k is the Boltzmann 
constant, Z is the mean ionic charge, h is the Planck con¬ 
stant and Ue and Ui are the electron and ion particle densi¬ 
ties respectively. The free-free emission SED is independent 
of the frequency at low energies and has an exponential cut¬ 
off at a characteristic frequency, determined by the temper¬ 
ature of the gas. This cutoff frequency shifts towards higher 
values as the temperature of the gas increases. The total flux 
on the other hand is a strong function of the density of the 
emitting gas parcel jl^ oc p°. Therefore, the temperature 
of the gas decides the hardness of the spectrum while its 
density controls the overall luminosity. 

In order to study the effect of the X-rays from clusters 
on gas cooling rates, it is important to match their observed 
luminosities. The relation between the X-ray luminosity and 
the mass of the host halo has been extensively studied all the 
way from galaxy groups (Bharadwaj et ar]|2014l to cluster 


scales ( Vikhlinin et al.||2009 Pratt et al.||2009 Planck Col 


laboration et al.|2013 Wang et al. [2014 1. More recently An¬ 

derson et al. (20141 analyzed a sample of ~ 250000 locally 


brightest central galaxies selected from the Sloan Digital Sky 
Survey (SDSS) and stacked the X-ray emission from these 
haloes, as a function of the stellar mass of the central galaxy, 
using data from the ROSAT All-Sky Survey. They constrain 
the Lx — M^, relation in a wide range of central galaxy stel¬ 
lar masses shown in Fig. (blue points). These values are 
quoted for hot gas emission between 0.15 — 1.0 Rsoo of the 
cluster. 
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Figure 3. X-ray (0.5-2.0 keV) luminosities of clusters. The blue 
points are observational measurements taken from |Anderson et al.| 
) |2014[ |. The red square {hl4) and triangle {hl5) are the lumi¬ 
nosities of the synthetic clusters calculated using Eqs. |6|7| These 
luminosity measurements are made between 0.15 — 1.0 Rsoo- 


We construct two idealised haloes containing gas and 
DM, one of mass (Mhaio) Mq (hl4) and one of mass 
10^® Mq (hl5). The DM density profile of the these haloes 
follows a |Navarro et al.|fi997| (NFW) distribution. The con¬ 
centration parameter is set by the c^ir — Mhaio relation given 


in Dutton & Maccio (20141. The gas fraction of the haloes 
is chosen to match observations (Giodini et al.pOOQ Gon- 


|zalez et al.|2013| [Sanderson et al.|2013[ ). The gas in the halo 

is placed in hydrostatic equilibrium with the DM and the 
temperature of the gas is set to the virial temperature of the 
halo. The properties of these haloes are shown in Table 
These haloes are then evolved, using AREPO, with cooling, 
gravity and hydrodynamics turned on for 200 Myrs. The lu¬ 
minosities are then calculated using the following relations 
(see also Eq. 2 of Grain et al.||20l6 1 : 


ju,a^a 

= Ao,a{Ta) 

= Ao,a{Ta) 


\ : 
P- ] 


Xe pa 

{Xe -b Xi)2 ^rup ^rup ’ 


Va 


( 6 ) 


where Va is the volume of the a*^ gas cell, Xe = Ue/nn, 
Xi = Ui/nn, Pa is the density of the cell, nip is the pro¬ 
ton mass and nia is the mass the gas cell. The electron and 
ion densities of any given gas cell are calculated assuming a 
primordial composition. The total X-ray luminosity is then 
given by 


Halo 

Mhaio [Mq] Ndm 

-^gas 

Cvir 

IDM [kpc] 

e9as[kpc] 

/gas 

hl5 

10^® 

50000 

100000 

5.21 

7.57 

2.92 

0.10 

hl4 

lO^-* 

50000 

100000 

6.56 

3.57 

1.27 

0.08 

h511 

5 X 10^1 

20000 

40000 

8.0 

0.82 

0.37 

0.15 

hl5hr 

10^® 

2.5 X 10® 

5 X 10® 

5.21 

1.30 

0.60 

0.10 

h2121r 2 X lO^^ 

5000 

10000 

10.0 

1.30 

0.60 

0.15 


Table 1. Table listing the parameters of the simulated galaxies 


comparison to the observational data points from | Anderson| 
et al. (20141 in Fig. The stellar mass of the haloes in 
question are determined by the M* — Mhaio relation given 
by Kravtsov et al. ( |2014[ ). The X-ray luminosity of the syn¬ 
thetically constructed clusters matches quite well the obser- 
vationally derived values. 

Using the relations given by Eq.[^ we then calculate the 
hot gas free-free emission fluxes at 100 kpc (red curve), 300 
kpc (green curve) and 1 Mpc (blue curve) from the center of 
the hl5 cluster (shown in Fig.|^. The SEDs show the char¬ 
acteristic shape of free-free emission, with more than one 
temperature component contributing to the overall shape. 
At 2 = 0, the hot halo cooling flux is much larger than the 
contribution from the extragalactic UV background (purple 
curve). 

The cooling curves of a gas parcel of nn = 0.001cm“® 
and Z = Zq under the influence of the hl5 cluster hot-halo 
radiation fields are shown in Fig. The reduction of cool¬ 
ing between 10"^ K and 10® K is quite drastic even at a dis¬ 
tance of 1 Mpc, with Fe being the only ionic species not 
completely photoionized. By comparison, a gas parcel at 
100 kpc is completely ionized. The only cooling contribution 
is from bremsstrahlung emission. The equilibrium tempera¬ 
ture, i.e., the temperature at which the heating and cooling 
rates of the gas parcel are equal, under these conditions also 
increases to about 3 x 10® K. 

Another important point to note is that the cooling 
rate of gas with temperatures greater than ~ 10®'® K does 
not change and seems independent of the incident radiation 
field. The high temperature plasma is completely thermally 
ionized and the only cooling mechanisms are bremsstrahlung 
and Gompton cooling. As the virial temperature of clusters 
are also above ~ lO^K, most of the ICM will be optically 
thin to the X-ray emission, which is in agreement with ob¬ 
servations ( Sarazin|1986 and references within). Therefore, 
the radiation field will not affect a large fraction of the clus¬ 
ter gas, but it will certainly reduce the cooling rate of any 
gas in the cluster with relatively low temperatures, either 
present in the form of cooling flows or as part of satellite 
galaxies. We can therefore conclude that the hot halo radi¬ 
ation fields will have a non-negligible effect on gas cooling 
rates in clusters. 


3 CALCULATING THE RADIATION FIELD 


^2 keV 

L = V / Lo,adn. (7) 

^ Jo.5 keV 

The luminosity calculated using these relations for the clus¬ 
ters hl4 (red square) and hl5 (red triangle) is shown in 


In the previous section the different sources of X-ray radi¬ 
ation were enumerated and the hot ICM cooling radiation 
was shown to be important. In RK14b we used an optically 
thin approximation, with a fixed escape fraction for low en¬ 
ergy photons. This allowed us to mimic absorption and at 
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Figure 4. The incident photon flux as a function of energy due to 
bremsstrahlung emission from the hot ICM, at various distances 
from the center of the hl5 cluster. The cyan curve is the flux at 
100 kpc using the three temperature approximation described in 
Section]^ The other curves show the fluxes at 100 kpc (red curve), 
300 kpc (green curve) and at 1 Mpc distance from the center of 
the cluster, calculated using the exact relations given by Eq. 
For comparison, the estimate for the flux of the extragalactic UV 
background at z = 0 is also plotted (purple curve). 


the same time make the radiative transfer problem compu¬ 
tationally feasible. 

The low energy part of the XRB spectrum is not ob¬ 
servable due to absorption by the disk of the host as well 
as our own Milky-Way (MW) galaxy. This necessitates that 
we only model XRB emission in the hard X-ray regime. The 
interaction cross section for the hard X-ray photons is small 
and hence their escape fraction from the disk of the galaxy 
is expected to be unity. Therefore the optically thin approx¬ 
imation holds for the XRB radiation fields. Additionally, we 
have argued in Section [2^ that most of the hot inter cluster 
gas is optically thin to the radiation field. This approxima¬ 
tion becomes invalid at the centers of the clusters, where the 
gas becomes dense enough to self-shield against the incident 
radiation field. This effect is modelled by assuming that gas 
above a certain density no longer sees the incident radiation 
held (see Sectionj^for more details). Therefore the optically 
thin approximation is valid for both radiation fields consid¬ 
ered in this paper making radiation transfer a simple inverse 
distance squared problem. 

3.1 Three temperature approximation 

The RK14b model, made a couple of simplifying assump¬ 
tions in order to quantify the effect of local radiation fields 
on gas cooling rates in galaxy formation simulations. One, 
except for absorption at the source, the gas is assumed to 
be optically thin at all wavelengths and two, the shape of 
each contributing component (young stars and post-AGB 
stars) to the overall SED does not vary with the underly¬ 
ing properties of the component. For example, the model 
assumes that the SED from old stars is constant in shape 



Figure 5. The volumetric cooling (solid curves) and heating (dot 
dashed curves) rates for a gas parcel under the influence of hot 
halo radiation field from the hl5 cluster at a distance of 100 kpc 
(red curves), 300 kpc (green curves) and 1 Mpc (blue curves). 
We also show the cooling and heating rates for an SED at 100 
kpc calculated using the three temperature approximation (cyan 
curves) described in Section For comparison the cooling and 
heating rates in the presence of an extragalactic UV background 
(purple curves) at z = 0 is also shown. The gas is at a density of 
nu = 0.001cm“®, and a metallicity of 7iq. 

irrespective of the age of the stellar population (as long as 
t* ^ 200 Myrs). This made sure that the spectrum of each 
contributing component does not change shape, rather only 
the overall normalisation of the spectrum varied according 
to the availability of and distance from the source. Therefore 
the photoionization (T) and photo heating (e) rates of each 
atomic species ‘i’ present in the gas, were just the sum of 
the individual components ‘c’, weighted according to their 
relative flux contributions ‘0’ 

(Ti, €i) = ^ (pcirc.i, Ec.i). (8) 

c 

This same method can be used for the XRB emission, be¬ 
cause they scale exactly in the same way as the young stars 
SED (which scales with SFR) and and emission from post- 
AGB stars (which scales with the total stellar mass of the 
galaxy) considered in RK14b. 

However, it is more difficult to model the emission from 
the hot gas. As discussed in Section [2.2[ the temperature of 
the emitting gas parcel controls the hardness of the spec¬ 
tra while the density controls the overall flux. The emitted 
spectra of a cluster will vary depending on the gas tempera¬ 
ture and density distribution within the cluster. Therefore, 
in addition to the overall normalisation, the shape of the 
SED also changes. The hot halo emission, therefore, cannot 
be quantihed as one of the components of Eq. instead one 
would theoretically need an infinite number of components, 
spanning the entire temperature range of the ICM. Every 
component of the hot halo radiation field will contribute an 
additional dimension to the cooling table. 

Observationally, the most dominant sources of hot gas 
emission are galaxy groups and clusters. In lower mass haloes 
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the hot halo emission is orders of magnitude lower and sub¬ 
dominant with respect to other X-ray sources like XRBs 


(Bogdan et al.|2013a|b Anderson et al.|2014 Bogdan et al. 


2015). Since galaxy groups/clusters usually have virial tem¬ 
peratures T > 10® ®K, we consider only cooling radiation 
from relatively hot gas; i.e. gas with logio(T[K]) > 5.5. 
We further sub divide this hot gas into three temperature 
regimes, 10®K, lO^K and 10®K, with a spread of 0.5 dex in 
each direction around each temperature bin. The gas parti¬ 
cles in each bin is then assumed to have the representative 
temperature of that bin. The flux from hot halo emission at 
any point ‘p’ is then approximated as 

= ’ ie [6,7.8], (9) 


(Xe-bW)2 47r^ 


where 




P atria 




V particles ‘a^. 


( 10 ) 


such that 


t-0.5<logio(Ta[K])^t-b0.5. (11) 

This approximation reduces the number of components 
of the hot halo emission radiation to three. The SED of each 
component has a fixed shape due to the hxed temperature 
and the normalisation just depends on the density times 
the mass of the particles within that temperature range. To 
check the validity of this approximation we plot, in Fig. 
the radiation flux at a point 100 kpc from the center of the 
hl5 cluster using the approximation (cyan curve) compared 
to the flux at the same point using exact relations (red curve) 
described in Section |2.2| The approximated SED matches 
quite well with the actual SED, with slight differences espe¬ 
cially in the regime where the SED changes shape, around 
100 Rydbergs. The corresponding cooling curve is plotted 
in Fig. The heating rates (dot dashed lines) lie on top of 
each other, while there seems to be a small difference to the 
cooling component. This difference is however of the order 
of 0.05 — 0.1 dex. Compared to the computational simplicity 
it offers, the errors induced by the approximation are quite 
small. 


3.2 Propagating the radiation field 

Now that we have found an effective way to quantify XRBs 
and hot halo emission, we need to come up with a compu¬ 
tationally efficient method to propagate the radiation held 
over the whole simulation volume. We have argued that the 
radiation helds considered in this work, for all practical pur¬ 
poses, can be assumed to be optically thin. This allows us 
to quantify the hux from each component in the parameter 
“(()’, which is proportional to the availability of the sources 
and inversely proportional to square of the distance between 
the source and the sink. 

We use the same method described in RK14b to cal¬ 
culate the hux of the radiation held; i.e. we tag on to the 
gravity tree of the simulation code. The tree structure par¬ 
titions the mass distribution into a hierarchy of localised 
regions. The force from nearby particles is calculated accu¬ 
rately, while the distant regions are explored more coarsely 
by treating them as single massive pseudo-particles, cen¬ 
tered on the center of mass of the corresponding region. 


In RK14b we implemented the LPF model in GASOLINE 
( Wadsley et al.|2004 | which used a k-D (binary) tree struc¬ 
ture ( Stadel||2001 1 in which bisections are done recursively 
through their longest axis of the cells, whereas, in this paper 


we implement the model in AREPO (Springel 20101 which 


uses the same oct-tree (Barnes & Hut 19861 structure as 
Gadget-2 ( Springel et al.|2005 |. Although these codes dif¬ 
fer in the actual construction of the gravity tree, they are 
the same for our purposes. 

One very simple improvement we have made, while tag¬ 
ging on to the gravity tree in Arepo is to consider an effec¬ 
tive center mass of the node while calculating the distance 
between the node and the sink position. The effective cen¬ 
ter of mass of a radiation field component ‘c’ in any node is 
defined as the center of effective mass of all the particles in 
that node that contribute to that particular component. The 
effective mass of a component in turn is defined as quantity 
that controls the normalisation of the SED of that com¬ 
ponent. For example, effective mass of the new stars SED 
component is just the mass of stars less than 10 Myr old (see 
RK14b for more details), while the effective mass of the hot 
halo emission components is the the mass times the density 
of the gas particles. This assures that the effective center of 
mass of a component in the node is close to the point in the 
node where the contribution to that particular component 
is maximum. 


4 COOLING TABLE CREATION 


We tabulate the gas cooling rates under various gas den¬ 
sities, temperatures and incident radiation fields. The to¬ 
tal gas cooling rate is divided into primordial, metal and 
Compton cooling (Eq. 15 of RK14b). The primordial and 
Compton cooling rates are calculated on-the-fly (see Vogels- 
berger et al.|2013] for more details). The metals are assumed 
to be in ionization equilibrium and their heating and cool¬ 
ing rates are tabulated across a range of physical conditions 
using CLOUDY. This look-up table is then used in simula¬ 
tions. 

Above 2 = 9, the gas is assumed to be of purely primor¬ 
dial composition, negating the need for a metal cooling table. 
From 2 = 9 to 2 = 3 the cooling table has five dimensions 
corresponding to the temperature, density, flux from new 
stars (t* ^ 10 Myrs), flux from old stars (t* ^ 200 Myrs) 
and redshift components. At 2 < 3 the table has eight di¬ 
mensions, the additional dimensions corresponding to the 
three hot halo emission components, which start to become 
dominant only at low redshifts once sufficiently massive 
haloes start forming. 

The density ranges from 10“® to 10^ cm“^ with a spac¬ 
ing of 1.0 dex in log space. The temperature ranges from 10^ 
to 10® K with a resolution of 0.1 dex. 4>sfr, which quanti- 
hes radiation from young stars and the HMXBs, varies from 
10“® to 10® Mq yr“® kpc“®. (jios, which quantifies the ra¬ 
diation field from post-AGB stars and LMXBs, varies from 
10® to 10®® Mq kpc“®. The spacing for both these 4>''s is 
1.0 dex. All three hot halo components, 4>t7 & 
ranges from 10®®"'® to 10®®'® cm“®, with a spacing of 2.0 dex. 
The redshift dimension that accounts for the change in UV 
background ranges from 0.0 to 9.0 with a resolution of 0.5. 

To create the table, cooling and heating rates are cal- 
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Figure 6. The temperature evolution of gas parcels with densities nn = 0.001 cm“^ (left panel) and nn = 0.1 cm“^ (right panel) and 
temperatures 3 X lO^K (solid curves), 3 X 10®K (dashed curves) and 3 X 10®K (dot dashed curves). Each gas parcel is subjected to two 
different radiation fields, one which contains hot halo radiation field of strength i(>t 7 = 10^^ cm“® plus the extragalactic UV background 
at 2 : = 0 (green curves) and other which contains only the extragalactic UV background flux at 2 : = 0 (blue curves). 


culated at every point of phase space for solar and primor¬ 
dial metallicity gas using CLOUDY. The difference between 
the solar and primordial metallicity values are stored as the 
heating and cooling rates due to metals only. 


5 TEST PARTICLE EVOLUTION 


As a first test of the cooling implementation in AREPO, 
we calculate the temperature evolution of isolated gas cells 
with and without hot halo radiation fields (Fig.|^. We con¬ 
sider gas cells of density of uh = 0.001 cm“® (left panel) 
and Uh = 0.1 cm“^ (right panel), metallicity Z = Zo and 
starting temperatures 3 x lO^K (solid curves), 3 x 10®K 
(dot dashed curves) and 3 x 10®K (dot dashed curves). The 
temperature evolution of these six gas cells under the the 
influence of a hot halo radiation field of strength (fiT 7 = 
10^^ cm“® plus the UV background at 2 ; = 0 (green curves) 
and under the influence of a UV background radiation field 
at z = 0 only (blue lines) are shown. 

The effect of the radiation field on the cooling rates 
of the gas cells is highly dependent on both the starting 
temperature and density of the cell. At a particular den¬ 
sity, the initial temperature of the gas decides the effec¬ 
tiveness of the radiation fields. The temperature evolution 
of the 3 X 10^ K gas cell (at both densities) is barely ef¬ 
fected. The time it takes for this high temperature gas cell 
to cool from 3 x lO^K to 10®K is much longer (~ 20 Gyrs for 
Uh = 0.001 cm“®) than to cool from 10®K to the equilibrium 
temperature (~ 1.7 Gyrs for uh = 0.001 cm“®). Therefore, 
even if the radiation fields drastically increa ses t he cooling 
time from 10®K to 10"^K (as seen in Section 2.21, the total 
cooling time will not change by a lot. 

Another notable change between the temperature evo¬ 


lution of gas cells with and without the hot halo radiation 
fields is the increase in the minimum temperature for each 
particle from ~ 10‘*K to ~ 10®K. The minimum tempera¬ 
ture reflects the equilibrium temperature the radiation field 
sets. While 10® K is a higher temperature than that at which 
stars form, it is not high enough to support against collapse 
with thermal pressure in cluster environments. As we will 
see in next section, 10®K gas collapses into a center that 
eventually reaches a high enough density to cool regardless 
of the radiation field. 

The cooling rates of the low temperature gas parcels 
(3 X 10®K & 3 X 10®K) on the other hand are reduced quite 
drastically. In the presence of the hot halo radiation field 
the equilibration timescale is about 7 times higher for 3 x 
10®K gas parcel and 25 times higher for 3 x 10®K gas parcel 
at Uh = 0.001 cm“®. In galaxies with < 10®'®K, the 
cooling and star formation timescales are inherently smaller. 
Therefore an increase in the cooling time of the gas particle 
will definitely affect the star formation rates of the galaxy. 
This implies that any galaxy that is close to or is falling into 
the cluster will experience substantial quenching, due to the 
host halo radiation field. 

The density of the gas also plays a major role in its 
cooling rates. As we discussed in Section [2.2[ the volumet¬ 
ric cooling rate is proportional to the square of the density, 
which in turn means that the cooling time of the gas par¬ 
cel will be inversely proportional to the density of the gas 
Ucooi = 3nkBT/2A). This affect can clearly be seen in Fig. 
M where the cooling times of gas cells with nn = 0.1 cm“®, 
irrespective of their starting temperature, is about a 100 
times shorter than the cells with nu = 0.001 cm“®. This 
means that galaxies with higher mean gas densities will have 
a larger cooling flow. 

Additionally the effect of the radiation fields will also 
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reduce at higher densities. We can define an ionization pa¬ 
rameter (U) for the gas parcel as 


U = 




( 12 ) 


where c is the speed of light, h is the Planck constant, n is the 
particle density and vt is the threshold frequency for ioniza¬ 
tion (13.6eV for hydrogen atom). This parameter is basically 
the ratio between the number density of ionizing photons in¬ 
cident on the gas cell to the number density ionic species in 
that cell. For a fixed radiation field the ionization parameter 
decreases as the density of the cell increases. A decrease in 
the ionization parameter, implies a decrease in the effective¬ 
ness of the incident radiation field. This effect is seen Fig.|^ 
where the relative difference in cooling times of the gas cells 
under the presence and absence of hot halo radiation fields 
is much smaller for the high density gas cells. For example, 
the gas cell with T^tart = 3 x lO^K (solid curves) has a three 
fold increase in the cooling time due to the hot halo radi¬ 
ation field when it has a density of nn = 0.001 cm~^, but 
the cooling times are practically identical when we increase 
the density of the cell to nn =0.1 cm“®. This means that 
galaxies with lower mean densities will be more amenable to 
reduction in cooling rates due to local radiation fields. 



logio{T[K]) 


Figure 7. A 2-D histogram of the amount of gas mass in the hl5 
cluster plotted as a function of the change in cooling rate between 
the UVB and UVB-I-LPF simulation runs and temperature at the 
end of 0.5 Gyr of evolution. 


6 RESULTS 

In order to study the effect of the hot halo radiation fields in 
a full dynamical setting, we perform isolated galaxy simnla- 
tions. Specifically, we simulate the idealised clusters hl4 and 
hl5 (see Section ! 


2.2 


for more details about the construction 


and properties of these clusters). 


6.1 Host Halo Quenching 

We simulate two clusters of halo mass 10^"^ M© (hl4) and 
10^® Mq (hl5). The only physical processes considered are 
gravity, hydrodynamics and gas cooling. A very simple pre¬ 
scription for self shielding is also included by imposin g the 
condition that all gas particles with nn Ss 0.1 cm“® (Cev- 
|erino et al.||2010l ) have zero local radiation flux. All the gas 
in the halo is assumed to have a metallicity of O.3Z0, a value 
chosen to match the observational constraints from nearby 
clusters ( |Sato et al.|2009[ ). 

Each cluster is simulated twice, once with gas cooling 
rates calculated in the presence of the UV background ra¬ 
diation field at z — 0 only (UVB) and once with the hot 
halo radiation field turned on in addition to the UV back¬ 
ground (UVB-I-LPF). The simulations are run for 0.5 Gyr. 
The percentage change in the amount of cold gas (T ^ 10"^K) 
formed between UVB and UVB-I-LPF simulation for the hl4 
and hl5 clusters is rather minimal (< 3%). This is because 
the virial temperature of clusters are high enough to fully 
thermally ionize the gas, which reduces the effectiveness of 
the incident radiation fields. This is seen more clearly in Fig. 

which shows the amount of gas mass as a function of tem¬ 
perature and the change in cooling rate between the UVB 
and UVB-I-LPF simulations for the hl5 cluster. Most of the 
gas in the halo lies in the region with T > 10^ ®K, where the 
change in cooling rate is pretty minimal. There is a > 10% 


(denoted by the solid horizontal black line in the Fig.[^ de¬ 
crease in the cooling rate of the UVB-I-LPF simulation only 
in the region where T < 10^'®K. However, the amount of 
gas in this region of the plot is only a small fraction of the 
total gas mass. There is in fact a wide region between lO'^K 
and 10®K, devoid of gas, because the cooling time in this re¬ 
gion is much lower than the dynamical time of the system. 
Below lO^K the radiation field does have a large effect, but 
gas at this temperature will enter the modified equation of 
state according to the ISM model used in our simulations 


(Springel & Hernquist 20031, meaning that it has already 


been earmarked for star formation. Therefore, changing the 
cooling of this gas will not reduce the star formation rate in 
our simulations. 

Another way to quantify the effect of local radiation 
fields is to look at the change in cooling timescales as a func¬ 


tion of the free-fall time of the system. Voit et al. (20141 have 


shown that clusters with on-going star formation and/or 
AGN activity have cooling times about the order of pre¬ 
cipitation time scales in clusters, which are empirically seen 
to occur when tff < tcooi < 20tfi ( Gaspari et al.|20ll| 20121, 
whereas, clusters with no star formation or AGN activity 
were shown to lie above this precipitation zone, meaning 
some mechanism, other than SF or AGN, is responsible for 
stabilising the cooling in these haloes. Fig. shows the gas 
cooling times as a function of radii for both the hl4 (blue 
curves) and hl5 (green curves) cluster simulations. The solid 
lines show the cooling times for the UVB simulation and the 
dot-dashed lines show the cooling times for the UVB-I-LPF 
simulation. These values are shown for the cluster after 0.5 
Gyr of evolution. In the hl5 cluster the UVB and UVB-I-LPF 
cooling times lie on top of one another, however, the hl4 
cluster does show a modest increase, a factor of 2, in the 
cooling times in the center. There are a couple of reasons 
why the radiation field is marginally more effective in the 
hl4 cluster. One, the virial temperature of the hl4 cluster 
is lower, meaning that the photons are more likely to inter- 
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Figure 8. Cooling time plotted as a function of radius in the hl4 
(blue curves), hl5 (green curves) and h511 (red curves) simula¬ 
tions. The dashed lines are lOtg, and the corresponding shaded 
regions are htg to 20tff. The solid lines represent the cooling rates 
in the UVB simulation while the dot-dashed lines are the cooling 
rates in the UVB-I-LPF simulations. The vertical lines are r = 2.5e 
above which the properties of the galaxies are well resolved. 


act with the gas. Two. the free fall times in the hl4 cluster 
are slightly longer, meaning a lower central density, which 
in turn decreases the density of the gas. This decreased gas 
density increases the cooling times and will also be more 
amenable to reduction in cooling rates due to local radia¬ 
tion fields (see Section |^. However, in neither cluster is the 
effect of the radiation field strong enough to increase the 
cooling timescales above the precipitation timescales of the 
system. 

We find that the low density gas present in the outskirts 
of the cluster has too high a temperature for the radiation 
field to effectively couple with it. On the other hand, the 
relatively cold gas, with large photon interaction cross sec¬ 
tions, are at too high a density, thereby reducing the effect 
of the radiation fields. We therefore conclude that hot halo 
radiation fields cannot significantly impact star formation in 
massive clusters. 


6.2 Satellite Quenching 

We have shown in Section|2^and Section[^ that the cooling 
rate of warm-hot gas (10"^ — 10®K) in clusters is significantly 
reduced in the presence of hot halo radiation fields. It follows 
that in systems where the dynamical timescales are similar 
to the cooling time of the warm-hot gas, the effect of the 
radiation fields will be much more pronounced. Low mass 
galaxies (Mhaio < IO^^Mq) have virial temperatures in the 
range of 10^ — 10®K and the dynamical timescales of the sys¬ 
tem are of the order of ~ 100 Myrs. This means that a hot 
halo radiation field, such as the one due to the ICM of the 
hl5 cluster, will change the dynamics of cooling flows in low 
mass systems. Therefore, satellite galaxies within the clus¬ 
ter will be adversely affected by this “radiation quenching" 
mechanism. 


Figure 9. Radiation field intensity plotted as a function of ra¬ 
dius for the hl4 (dot-dashed curves) and hl5 (solid curves) clus¬ 
ters separated into the three temperature components tj)T6 (blue 
curves), (/)x 7 (green curves) & 0x8 (red curves). 


To test this scenario, we construct a galaxy of mass 
5 X 10^^ Mq, labelled h511, in the same way as described 
in Section |2.2| The properties of this galaxy are given in 
Table The virial temperature of this halo is 10® ®K. 
As shown in Section |2.2| the radiation field of the cluster 
varies with the distance from the center (Fig. |^. There¬ 
fore the position of the satellite within the cluster will de¬ 
termine how effective the host radiation fields can be. Fig. 

shows flux in each component of the hot halo radiation 
field, (0x6 blue curve, 0 x 7 green curve, 0x8 red curve), as 
a function of distance from the center of the corresponding 
cluster (hl4 - dot dashed curves and hl5 - solid curves). The 
virial temperature of the hl5 cluster is ~ 10®K, therefore the 
amount of gas in the 10®'® < T ^ 10®'®K range is negligible 
and hence the 0x6 component is sub-dominant. Whereas, 
the virial temperature of the hl4 cluster is ~ lO^K, hence 
the 0x8 component is non-existent. 

We place the h511 halo at a distance of 100 kpc (blue 
curves) and 300 kpc (green curves) in both the hl4 (dot 
dashed curves) and hl5 (solid curves) clusters and simulate 
the galaxy for 0.5 Gyrs with the host halo radiation field 
turned on. 

There is ~ 40% reduc tion in the amount of cool 
(T < 10"^K) gas formed (Fig. 101, if the satellite is placed 
at a distance of 100 kpc from the center of the hl5 cluster 
compared to the field and the effect reduces to ~ 20% at the 
distance of 300 kpc from the center. The effect is smaller if 
the galaxy is placed in the hl4 halo with the reduction in 
amount of cool gas about 10%. 

Now the question arises as to why we see an affect on 
low mass galaxies. To answer this we plot the gas mass of 
the satellite as a function of the fractional change in the 
cooling rate of gas between the UVB and UVB-bLPF sim¬ 
ulations and its temperature, when placed at a distance of 
100 kpc from the center of the hl5 cluster (Fig. 111. The 
cooling rates of almost all the gas in the satellite is severely 
reduced in the presence of the host radiation field. This in 
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Figure 10. Percentage change in the amount of cold gas (T < 
lO^K) formed between the UVB+LPF and UVB simulations as 
a function of time in the h511 galaxy simulation when placed at 
100 kpc (blue curves) and 300 kpc (green curves) from the center 
of the hl5 (solid curves) and hl4 (dot dashed curves) clusters. 



^_____ i ■ 

4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 


logio{T[K]) 

Figure 11. A 2-D histogram of the amount of gas mass in the 
h511 halo plotted as a function of the change in cooling rate be¬ 
tween the UVB and UVB-I-LPF simulation runs and temperature 
at the end of 0.5 Gyr of evolution. The h511 galaxy is assumed 
to be at a distance of 100 kpc from the center of the hl5 cluster. 


turn causes a reduction in the amount of cooling flow and 
hence a reduction in the amount of cold gas formed. In con¬ 
trast the amount of gas with reduced cooling rates due to 
the hot halo radiation field was minimal for the host cluster 


(see Fig. 
Fig. 


shows the cooling time of the gas in h511 halo, 
at a distance of 100 kpc from the center of the hl5 clnster, 
as a function of radii (red curves). In the absence of local 
radiation fields the cooling time scales (red, solid curve) are 
~ 10“"* tprep (red, dashed curve and shaded region), imply¬ 
ing a high star formation rate. The radiation fields increase 
the cooling times by a factor of about 100 — 1000 (red, dot- 
dashed curve), making them only a order of magnitude lower 
than the precipitation timescales. This huge effect is due to 
the fact that, unlike clusters, most of the gas in the satel¬ 
lite is at exactly the temperatures at which the ionization 
cross section in maximum (< 10®K), and at low densities, 
as evidenced by the constantly higher free fall times in the 
satellite. These increased cooling times will lead to a drastic 
reduction in star formation rate of the galaxy, just because 
of lower amounts of cooling flows (see RK14b). 


6.2.1 Relative effectiveness of the radiation quenching 
mechanism 

Radiation quenching complements other environmental 


quenching mechanisms such as ram pressure ( 

Gunn & Gott 

19721 IMcCarthy et al.||20081 |Font et al.||2008 

|Kang & van 

den Bosch||2008 

Weinmann et al.||2010 

1 and tidal stripping 

(Taylor & Babul 20011. Therefore, it 

is important to un- 


derstand the impact of radiation quenching in the presence 
of these other environmental quenching mechanisms. |Mc-| 
Carthy et al. (20081 have found that for satellite galaxies 


with typical structural and orbital parameters, up to 30% of 
the initial hot halo gas can remain with the galaxy for up to 
10 Gyr. During this phase it was assumed that the halo gas 


will continually cool and replenish the disk and sustain star 
formation in the satellite. However, our results show that 
the host halo radiation fields will reduce the star formation 
in these satellites to a certain extent, even before they are 
completely stripped of gas. 

To test this scenario we run a non-cosmological merger 
simulation of a 2 x 10^^ Mq halo falling into a 10^® Mq 
cluster. We construct a high resolution version of the hl5 
cluster, where the DM distribution is sampled by 2.5 mil¬ 
lion particles and the gas distribution sampled by 5 million 
particles (hl5hr). The galaxy of mass 2x 10^^ M© (h2121r) is 
sampled with 5000 DM particles and 10000 gas cells. These 
galaxies are constructed in the same way as described in 
Section [2.2| and the properties of them are listed in Table 
The h2121r satellite galaxy is placed at the virial radius of 
the hl5hr cluster with the velocity of the satellite given by 
the relations outlined in Benson ( 2005| ). The system is then 
evolved for 1 Gyr with (UVB-I-LPF) and without (UVB) the 
hot halo radiation fields turned on. In addition to gravity, 
hydrodynamics and gas cooling, we also turn on star forma¬ 
tion and stellar feedback according to the models described 


Vogelsberger et al. (20131. 


Fig. 12 shows the amount of stellar mass (solid curves) 


formed and the instantaneous star formation rates (dashed 
curves) in the h2121r galaxy as it falls in and orbits around 
the hl5hr cluster. The amount of stellar mass formed af¬ 
ter 1 Gyr is about 20% lower in the UVB-I-LPF run (red 
curve) compared to the UVB run (green curve). Addition¬ 
ally, the star formation rate in the UVB-I-LPF simulation is 
also about 20—30% lower. These results suggest that even in 
the presence of other environmental quenching mechanisms, 
radiation quenching plays a significant role in reducing star 
formation and stellar mass in cluster environments. There¬ 
fore, they can in principle account for the reduced star for¬ 


mation rates of galaxies in cluster environments (Poggianti 
et al. 1999 Mei et al.|2009 1 and also explain why local qui- 
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Figure 12. The amount of stellar mass formed (solid curves) 
and the instantaneous star formation rate (dashed curves) in the 
UVB+LPF (red curves) and UVB (green curves) runs as a func¬ 
tion of simulation time in the h2121r satellite halo, when orbiting 
the hl5hr cluster. 


escent galaxies are observed to be quenched by a ” strangu¬ 
lation “ mechanism rather than by sudden outflows of gas 
due to SNe feedback or ram pressure stripping (|Peng et al.| 
201^ |Lu et 


7 CONCLUSIONS 

In this paper, we present first simulations of galaxy forma¬ 
tion which include the impact X-ray emission from binaries 
and hot ICM on gas cooling rates using the AREPO code. 

We find the observed X-ray flux from XRBs to be sub¬ 
dominant compared to other stellar radiation sources, min¬ 
imally affecting the gas cooling rates in galaxies. However, 
hot gas emission from clusters has a high photon flux all the 
way to the virial radius of the halo. This is because the ICM 
is an extended source, therefore the bremsstrahlung emis¬ 
sion from the ICM will not fall off as rapidly as radiation 
from centralised sources like stars. The gas cooling rates are 
substantially reduced, especially in the lO^K < T < 10®K 
regime, under the influence of these radiation helds. 

We test our model in a full dynamical setting by sim¬ 
ulating idealised clusters of mass 10^'* Mq and 10^® M©. 
The difference between the amount of cold gas (T < lO'^K) 
formed between the simulations with and without local ra¬ 
diation feedback is minimal, because the cooling flows are 
not affected by the local radiation sources in clusters. This 
is because the low density gas present in the outskirts of the 
cluster is at too high a temperature for the radiation field to 
effectively couple with it. On the other hand, the relatively 
cold gas, with large photon interaction cross sections, are at 
too high a density, which reduces the effectiveness of the ra¬ 
diation fields. This means that, even though there is enough 
energy in the form of radiation to quench gas cooling in clus¬ 
ters, the coupling between the gas and the radiation field is 
too weak. We conclude that the radiation fields considered 
in this paper will not have a major effect on star formation 
rates in massive galaxies. 


The cooling rate of relatively cold gas (T < 10®K) 
on the other hand is reduced considerably by the hot-halo 
radiation field. The change in the cooling rate can be as 
large as the 30 times the fiducial value. This can in prin¬ 
cipal change the dynamics of cooling flows in the systems 
where the dynamical timescales are comparable to the cool¬ 
ing timescales of gas. This implies that relatively low mass 
satellites (Mhaio lO^^M©) will experience quenching due 
to the radiation field of the host halo in cluster environ¬ 
ments. We tested this idea by simulating a low mass galaxy 
(Mhaio = 5 X 10^^ Mq) in the presence of host halo radiation 
helds. We placed the satellite galaxy at a various distances 
from the center of the cluster and assumed that all gas in 
the satellite sees a host halo radiation held corresponding to 
that distance. These simulations show an apparent reduc¬ 
tion in the amount of cold gas formed when the host halo 
radiation held is turned on. At a distance of 100 kpc from 
the center of the halo, the satellite has a ~ 40% reduction in 
cold gas mass in 0.5 Gyrs. This reduces to about 20% at 300 
kpc and reduces further at larger radii. This effect is due to 
the fact that, unlike clusters, most of the gas in the satel¬ 
lite is at exactly the temperatures at which the ionization 
cross section is largest (T < 10®K), and have low enough 
densities, to be highly inhuenced by the radiation held. 

We also compared the ehectiveness of radiation quench¬ 
ing to other environmental quenching mechanisms (e.g., ram 
pressure stripping) by running a galaxy merger simulation of 
a 2x 10^^ Mq halo falling into 10^® M© halo. The UVB-I-LPF 
simulation shows a 20% reduction in the stellar mass of the 
satellite after 1 Gyr with respect to the UVB-only simula¬ 
tion. We also hnd a reduction of ~ 20 — 30% in star forma¬ 
tion rates at late times. These results indicate that even in 
the presence of other environmental quenching mechanisms, 
the hot halo radiation helds still play a signihcant role in 
quenching galaxies in cluster environments. 

We note that these results are preliminary and we plan 
to extend this initial study to a full cosmological simulation 
and quantify the ehectiveness of the local radiation helds 
on satellite quenching in a wide range of satellite and host 
galaxy masses. 
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